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Abstract 

Based on first-principles calculation and effective model analysis, we propose that the cubic 
antiperovskite material CusPdN can host a three-dimensional (3D) topological nodal line semimetal 
state when spin-orbit coupling (SOC) is ignored, which is protected by coexistence of time-reversal 
and inversion symmetry. There are three nodal line circles in total due to the cubic symmetry. 
“Drumhead”-like surface flat bands are also derived. When SOC is included, each nodal line evolves 
into a pair of stable 3D Dirac points as protected by C4 crystal symmetry. This is remarkably 
distinguished from the Dirac semimetals known so far, such as NasBi and Cd3As2, both having 
only one pair of Dirac points. Once C4 symmetry is broken, the Dirac points are gapped and the 
system becomes a strong topological insulator with (1;111) Z2 indices. 
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Introduction — Band topology in condensed matters has attracted broad interests in recent 
years. It was prospered by the discovery of two-dimensional (2D) and three-dimensional 
(3D) topological insulators (TIs) [H |2] • These materials exhibit a bulk energy gap between 
the valence and conduction bands, similarly to normal insulators, but with gapless bound¬ 
ary states that are protected by the band topology of bulk states. Topological classihcation 
has also been proposed for 3D metals [3j. The topological invariant, so called Fermi Chern 
number, can be dehned on a closed 2D manifold, such as the Fermi surface, in 3D momen¬ 
tum space [3HS]- This is essentially the same as Chern number defined in the whole 2D 
Brillouin zone (BZ) for insulators. Up to now, three types of nontrivial topological metals 
have been proposed. They are Weyl semimetal M, Dirac semimetal [U HO] and nodal line 
semimetal (NFS) [II11I2]. All of them have band crossing points due to band inversion [13] . 
For Weyl and Dirac semimetals, crossing points locate at different positions which com¬ 
pose Fermi surfaces. For the NFS, the crossing points around the Fermi level form a closed 
loop. The breakthrough in topological semimetal research happened in material realiza¬ 
tion of Dirac semimetal state in NasBi and Cd 3 As 2 , which were first predicted theoretically 
[H Ho] and then confirmed by several experiments [HHIZI. Starting from Dirac semimetal, 
one can obtain Weyl semimetal by breaking either time-reversal [ZIIHIIIHI or inversion sym¬ 
metry [T 9 H 22 ]. Among them, the prediction of Weyl semimetal state in nonmagnetic and 
noncentrosymmetric TaAs family [23l [2l| have been verihed by experiments [251130] . 

Existence of NFS has been proposed in Bernal graphite [HU [H2] and other all-carbon 
allotropes, including hyper-honeycomb lattices [33] and Mackay-Terrones crystal (MTC) [T^ . 
In Ref. dll it was discussed that when spin-orbit coupling (SOC) is neglected and band 
inversion happens, the coexistence of time-reversal and inversion symmetry protects NFS 
in 3D momentum space. It is heuristic that compounds with light element can also exhibit 
nontrivial topology, which is totally different from the common wisdom that strong SOC in 
heavy-element compounds is crucial for topological quantum states. The NFS state was also 
proposed in photonic systems with time-reversal and inversion symmetry [31]. It was also 
shown that mirror symmetry, instead of inversion symmetry, together with time-reversal 
symmetry protects NFS when SOC is neglected [33 [36] • TaAs [23], Ca 3 P 2 [3Z] and FaX 
(X=N, P, As, Sb, Bi) [3H] belong to this type. 

In the present work, based on hrst-principles calculations and effective model analy¬ 
sis, we demonstrate that the antiperovskite Cu 3 PdN is a new candidate for realizing NFS 
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and “drumhead”-like surface flat bands [T^ I391 - H3] . which may open an important route 
to achieving high-temperature superconductivity [HHIU]. Strong SOC will drive NLS in 
CuaPdN into Dirac semimetal state with three pairs of Dirac points, leading to exotic sur¬ 
face Fermi arcs which can be observed on various surfaces of this material. This is very 
unique since Dirac semimetals known so far, NasBi and Cd 3 As 2 , have only one pair of Dirac 
points. When the C 4 crystal symmetry in CusPdN is broken, the Dirac points are gapped 
and the system becomes a strong TI with Z 2 indices It is well known that the coop¬ 

erative interactions among lattice, charge, and spin degrees of freedom make antiperovskites 
exhibit a wide range of interesting physical properties, such as superconductivity mi, gi¬ 
ant magnetoresistance [IH], negative thermal expansion |l9j, and magnetocaloric effect |50] . 
Our prediction of 3D NLS and Dirac semimetal state in antiperovskites provides a promising 
platform for manipulating these exotic properties in presence of nontrivial topology. 

Crystal structure and methodology — Perovskite has a formula ABX 3 , where A and B are 
cations and X is an anion. The antiperovskite parallels to perovskite but switching the posi¬ 
tion of anion and cation, namely in ABX 3 X stands for an electro-positive cation, while A for 
an anion as shown in Fig. [^a) for CusPdN. Here we perform density functional calculations 
by using the Vienna ab initio simulation package [SI] with generalized gradient approxima¬ 
tion [S2] and the projector augmented-wave method [S3|. The surface band structures are 
calculated in a tight-binding scheme based on the maximally localized Wannier functions 
(MLWF) [Ml, which are projected from the bulk Bloch wave functions. 

Electronic structures — The band structure of Cu 3 PdN is shown in Fig. The fatted 
bands suggest that the valence and conduction bands are dominated by Pd-dd (blue) and 
Pd-5p (red) states. They also indicate band inversion at R point, where the Pd-5p is lower 
than Pd-4d by about 1.5 eV. To overcome possible overestimation of band inversion [55], we 
employ hybrid density functionals [56] to confirm the existence of band inversion, while the 
energy gap at R point is slightly reduced to 1.1 eV. 

Without SOC the occupied and unoccupied low energy bands are triply degenerate at R 
point (six-fold degenerate if spin is considered). These states belong to the 3D irreducible 
representations P 4 and P^ of Oh group at R point, respectively. We emphasize that, unlike 
the situation in typical TI such as the Bi 2 Se 3 family compounds EH, the band inversion in 
Cu 3 PdN is not due to SOC. To illustrate the band inversion process in Cu 3 PdN explicitly. 
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Figure 1; (Color online) (a) Crystal structure of antiperovskite CusPdN with Pm3m (No.221) 
symmetry. Nitrogen atom is at the center of the cube and is surrounded by octahedral Cu atoms. 
Pd is located at the corner of the cube, (b) Bulk and projected (001) surface Brillouin zone. The 
three nodal line rings (orange color) and three pairs of Dirac points (red points) without and with 
SOC included, respectively, are schematically shown. 

we calculated energy levels of and Fj bands at R point for CusPdN under different 
hydrostatic strains. As presented in Fig. |^b), the band inversion happens at a = l.lloo 
and the inversion energy increases as further compressing the lattice. The intriguing point 
of the CusPdN band structure without SOC is that the band crossings due to the band 
inversion form a nodal line circle because of the coexistence of time reversal and inversion 
symmetry as addressed in Ref. na. 

The protection of a nodal line can be inferred from the topological number of the form [3^ 

US] 

7 = y -A(k) ■ dk, (1) 

where ^(k) is the Berry connection of the occupied states, C is a closed loop in the momen¬ 
tum space. If C is pierced by the nodal line, one has 7=7r, otherwise 7=0. Below we prove 
the existence of nodal line from the argument of codimension. A general 2x2 Hamiltonian 
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Figure 2: (Color online) (a) Electronic band structure without SOC, where the component of Pd-Sp 
(4(i) orbitals is proportional to the width of the red (blue) curves. Band inversion between p and 
d orbits happens at R point, (b) Evolution of energy levels of E^ and EJ at R under hydrostatic 
pressure. Band inversion happens when a < l.llao. (c) Electronic band structure with SOC 
included. A small gap is opened in the R-X direction while the Dirac point in the R-M direction 
is stable and protected by crystal symmetry C 4 rotation. 

is enough for describing the two crossing bands, 

H(k) = gcc{k)acc + gy{k)ay + gz{k)a^, (2) 

where k ={kx, ky, kz), gx,y,zO^) are real functions and Cx^y^z are Pauli matrices characterizing 
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the space of the two crossing bands, which are mainly from and dxy orbitals of Pd in 
CuaPdN. The coexistence of time-reversal and inversion symmetry leads to Qy = 0, and 
Qz being odd and even fnnction of k [12]. Up to the lowest order of k, and gz(k) are 

given as 


gx(k) = jkxkykz, 

(3) 

gz{k)=M-B{kl + kl + kl), 

(4) 


considering the crystal symmetry at R point. The eigenvalnes of Eq. Q are U(k) = 
±'\/( 7 ^(k) -|- gliy). Nodal lines appear when 5 'a;(k) = 0 and gziy) = 0. It is easy to check 
that (y'^(k) = 0 can be satished only in the case of MB > 0, which is nothing bnt the band 
inversion condition. When band inversion happens, there alway exists three closed nodal 
lines in momentnm space as shown in Fig. [^b) which are the solntions of gx(^) = 5 'z(k) = 0. 

When SOC is considered, the six-fold degenerate states at R are split into one fonr-fold 
and one two-fold degenerate states. As shown in Fig. |^c), there are two fonr-fold degenerate 
states close to Fermi energy: the occnpied one with P^ symmetry and the nnoccupied one 
with Pg symmetry. Moving from R to X, the symmetry is lowered to C 2 v The fonr-fold 
degenerate states at R are split into two-fold degenerate states along R-X direction. The 
hrst-principles calcnlation shows that two sets of bands with the same Ps symmetry are 
close to the Fermi energy and a gap ~0.062 eV is opened at the intersection as shown 
in the inset of Fig. [^b). In order to rednce this SOC splitting and thus achieving NTS, 
one can replace Pd with lighter elements, such as Ag and Ni. Along R-M direction, the 
symmetry is characterized by C^y double group. As indicated in Fig. |^c), the two sets 
of bands close to Fermi energy belong to Py and Pg representation, respectively. They are 
decoupled and the crossing point on R-M path is unaffected by SOC. They form a Dirac node 
near the Fermi energy as shown in Fig. |^c), which is protected by the crystal symmetry 
C 4 rotation [a [ini EE]. Once C 4 rotational symmetry is broken, the Dirac node will be 
gapped and the parities of occupied states at eight time reversal invariant momenta (TRIM) 
point are shown in Table. The Z 2 indices are indicating a strong TI. The band 

structure of CuaPdN is different from that of antiperovskite SrgPbO (59] |60], where the 
band inversion happens at P point and the involved bands belong to the same irreducible 
presentation, which leads to anti-crossing along P-X direction. 

Effective Hamiltonian for the 3D Dirac fermions — In order to better understand the band 
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TRIM points 

R 

r 

M(x3) 

X(x3) 

parity product 

± 

- 

- 

- 


Table I: Parity product of occupied states at the TRIM points. The Z 2 indices are (1;111). 


crossing and gap opening discussed above, we derive an low energy effective Hamiltonian 
based on the theory of invariants in a similar way as that for 612803 family 

We hrst construct the model Hamiltonian along R-M in direction. The symmetrical 
operations in this direction contain the crystalline symmetry and inversion with time- 
reversal symmetry (PT). As discussed above, the wave functions of low-energy states along 
R-M direction are Pe symmetry with angular momentum = ±1/2 and Py with = ±3/2. 
Therefore, the model Hamiltonian respecting both and PT symmetries can be written 
as 


Hrm — 


A4.1 0 C 2 /C+ c^k'^ ± Cik"^ 

M.I C4k‘^_ + C'ik\ —C2k- 

M 2 0 

t M2 


(5) 


up to the second order of k in basis of {\j^ = \)p, \j^ = -\)p,\jz = \)d, \jz = “Dd}, where 
M-i = nip + milkz + mi 2 k‘^, M 2 = ±^ 21^2 ±^ 22 ^^, = kl + ky + k"^ and k± = k^^iky. 

The (74 rotation symmetry requires the matrix element between \jz = and \jz = — |) to 
take the form of /c^ in order to conserve the total angular moment along ^ direction. Due to 
PT symmetry and mirror symmetry, all parameters in Eq. ([^ can be chosen as real. Their 
values can be derived by htting the dispersions to those of hrst-principles calculations |63] . 

On kz axis, the effective Hamiltonian is diagonal and the | ± 1/2) sets and | ± 3/2) sets 
are decoupled. Since the energy of p and d orbitals are inverted at R point, namely mp < 
and 77122 < 0 < 77112, the | ± 1/2) and | ± 3/2) bands cross at = [((ttih — 77121 )^ — 4(77112 — 
m 22 ){mp - 777 rf))P 2 _ [rnii - m 2 i)\l[ 2 {mi 2 - 77722 )] point. 

The model Hamiltonian along R-X direction can be derived in the same way with C 2 V 
and PT symmetries. From the hrst-principles calculations, the states near Fermi energy are 
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characterized by Fs with = | ± 1/2). The model Hamiltonian can be written as 


Hnxik) 


Ml 0 AI3 cik- 
Mi cik+ —M3 
M 2 0 

t M2 


( 6 ) 


up to the second order of k in basis set of {\j^ = l>p, = -l>p, = l>s, = -l>s}. 

Here, k^ is set as along R-X direction. The term M 3 = m 3 + m 3 ikz + m 32 k‘^ conserves 
the total angular momentum along 2 ; direction. It couples the going-down Ts states and 
going-up Fs states and opens a gap at the crossing point shown in Fig. |^c). 


Surface states — The band inversion and the 3D Dirac cones in CuaPdN suggest the presence 
of topologically nontrivial surface states. They are calculated based on TB Hamiltonian 
from MLWF [5l IH]. The obtained band structures and surface density of states (DOS) on 
semi-inhnite (001) surface are presented in Fig. 

Without SOC, the bulk state is the same as MTC |T2] and there exists surface flat bands 
nestled inside the projected nodal line ring on the (001) surface, namely the “drumhead” 
states as shown in Figure|^a). The peak-like DOS from these nearly flat bands is also clearly 
shown. The small dispersion of this “drumhead” state comes from the fact that the nodal 
line ring is not necessarily on the same energy level due to the particle-hole asymmetry m 
[T 21 06] . This result is consistent with the correspondence between Dirac line in bulk and flat 
band at boundary as established in Ref. [391011 03] . Such 2D flat bands and nearly inhnite 
DOS are proposed as a route to achieving high-temperature superconductivity 0^106] . 

In the presence of SOC, the nodal line ring will be gapped in general. However, there 
is an exception. For example, in TaAs family each ring evolves into three pairs of Weyl 
nodes [23]. In CusPdN, each ring is driven into one pair of Dirac nodes. The (001) surface 
state band structure in Fig. |^b) clearly shows the gapped bulk state along f-M direction 
and the existence of surface Dirac cone due to topologically nontrivial Z 2 indices as seen in 
Na 3 Bi [4j and Cd 3 As 2 ra. The bulk band structure along R-X and R-M overlap each other 
when projected onto (001) surface along the X-M path. The bulk Dirac cones are hidden 
by other bulk bands. Therefore, it is difficult to identify the detailed connection of Fermi 
arcs in the Fermi surface plotting as shown in Fig. [^ though some eyebrow-like Fermi arcs 
can be clearly seen around these projected Dirac nodes. 




1 


(a) 


0.5 


> 

0) 

>• 

DJ 

0 ) 

c 


-0.5 


-1 


-1.5 

1 


0.5 



>. 

D) 

I -0.5 

UJ 


-1 


-1.5 




Figure 3: (Color online) Band structures and DOS for (001) surface (a) without and (b) with SOC 
. Without SOC, the nearly flat surface bands are clearly shown in the white dashed box around M 
point. 


Conclusion — In summary, we propose that 3D topological nodal line semimetal states can 
be obtained in a nonmagnetic and centrosymmetric system CusPdN. The “drumhead”-like 
surface flat bands nestled in a projected nodal line ring have been obtained. Including SOC 
will drive each nodal line ring into one pair of Dirac points to host Dirac semimetal state. 
The surface Dirac cone and the Fermi arcs around projected Dirac cones are observed. The 
existence of multiple pairs of 3D Dirac points distinguishes this system from other Dirac 
semimetals with only one pair. The cubic antiperovskite structure of CusPdN makes it 
a good platform for manipulating ferromagnetism, ferroelectricity and superconductivity 
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Figure 4: (Color online) Fermi surface of (001) surface shown in Fig. [^b) with chemical potential 
at 0.12 eV (green dashed line in Fig. [^(b). The yellow dots are projected Dirac points. The bigger 
one at M means there are two Dirac points superposed there. 

realized in a broad class of materials with perovskite structure in presence of nontrivial 
topology. 
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